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The collective variable representation (CV) of classical statistical systems such 
as simple liquids has been intensively developed by the Ukrainian school after 
seminal works by Prof. Ihor Yukhnovskii. The basis and the structure of the CV 
representation are reexamined here from the point of view of statistical field 
theory and compared with another exact statistical field representation of liquids 
based upon a Hubbard-Stratonovich transform. We derive a two-loop expansion 
for the grand potential and free energy of a simple fluid in both version of the 
theory. The results obtained by the two approaches are shown to coincide at 
each order of the loop expansion. The one-loop results are identical to those 
obtained in the framework of the random phase approximation of the theory of 
liquids. However, at the second-loop level, new expressions for the pressure and 
the free energy are obtained, yielding a new type of approximation. 
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1. Introduction 

Functional methods in modern statistical physics represent one of the most pow- 
erful tools for the study of both equilibrium and dynamical properties (see, e.g. [ 
QlEj). In many cases the partition function of simple models (Ising and Heisenberg 
spins or classical particles in interaction) can be re-expressed as a functional integral 
after performing a Hubbard-Stratonovich transformation, a simple device proposed 
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in the 50ies [ EH H] • Nearly at the same time another method - the so-called collec- 
tive variables (CV) method - that allows in a explicit way to construct a functional 
representation for many-particle interacting systems was developed [ |3J E] and ap- 
plied for the description of charged particle systems, in particular, to the calculation 
of the configurational integral of the Coulomb systems. The idea of this method is 
based on: (i) the concept of collective coordinates being appropriate for the physics 
of system considered (see, for instance, [Ej), an d (ii) an integral giving an exact 
functional representation for the configurational Boltzmann factor. Later the CV 
methods was successfully developed for the description of classical many-particle 
systems [|S] and for the theory of second order phase transitions The method 
has a wide range of application and, to paraphrase Pr. Yukhnovskii in his book [ 
|U] " The collective variables is a common name of a class of variables which are 
specific for each kind of physical system. For magnetic systems CV are variables 
corresponding to the modes of vibrations of spin moments, for ferroelectrics they 
are cluster vibrational modes, for systems of charged particles they are associated 
with the modes of generalized charge vibrations, for binary alloys they are modes 
of one-particle distribution function, the CV of the liquid-gas system are modes of 
some deviations from the critical density and so on." 

One of the goals of this paper is to reconsider the CV method from the point of 
view of statistical field theory and to compare the results obtained with those found 
recently by one of us by means of the KSSHE (Kac-Siegert-Stratonovich-Hubbard- 
Edwards) theory [ITU]. 

We formulate the method of CV in real space and consider a one-component 
continuous model consisting of hard spheres interacting through additive pair po- 
tentials. The expression for the functional of the grand partition function is derived 
and the CV action that depends upon two scalar fields - field p connected to the 
number density of particles and field u> conjugate to p and connected to the den- 
sity of energy - is calculated. We study correlations between these fields as well as 
their relations to the density and energy correlations of the fluid. The grand par- 
tition function of the model is evaluated in a systematic way using a well-known 
method of statistical field theory, namely the so-called loop expansion. It consists in 
expanding functionally the action TC around a saddle point, so that the lowest order 
(zero loop) approximation defines the mean-field (MF) level of the theory and the 
first order loop expressions correspond to the random phase approximation (RPA). 
Recently [EI] this technique was applied to the action obtained within the frame- 
work of the KSSHE theory. In this paper we perform a two-loop expansion of the 
pressure and the free energy of the homogeneous fluid which yields a new type of 
approximation which we plan to test in our future work. 

The paper is organized as follows. In Section (J2J), starting from the Hamilto- 
nian, we introduce the two different functional representations of the grand partition 
function based on the KSSHE and CV methods. Here we also enter several types 
of statistical field averages that are important in the further part of the paper. In 
Section (j3j) we introduce the CV and KSSHE field correlation functions, establish 
links between them as well as their relation to the density correlation functions of 
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the fluid. The MF level of the KSSHE and CV field theories is formulated in Sec- 
tion (HJ). Section (jSJ) is devoted to the loop expansion of the grand potential. The 
pressure and the free energy of the homogeneous fluid are obtained in the two-loop 
approximation in Section 6. We conclude with some final remarks in Section ©. 



2. Prolegomena 
2.1. The model 

We consider the case of a simple three dimensional fluid made of identical hard 
spheres of diameter a with additional isotropic pair interactions v (ry) (r^ = Xj\, 
Xi is the position of particle "z"). Since v (r) is an arbitrary function of r in the core, 
i.e. for r < a, we will assume that v (r) has been regularized in such a way that its 
Fourier transform v q is a well behaved function of q and that v (0) is a finite quantity. 
We denote by Q the domain of volume V occupied by the molecules of the fluid. 
The fluid is at equilibrium in the grand canonical (GC) ensemble, f3 = 1/k^T is 
the inverse temperature (k-g Boltzmann constant), and \x the chemical potential. In 
addition the particles are subject to an external potential if>(x) and we will denote 
by u(x) = (3 (fi — if) (x)) the dimensionless local chemical potential. We will stick to 
notations usually adopted in standard textbooks devoted to the theory of liquids 
(see e.g. [E]) and will denote by w(r) = —f3v(r) minus the dimensionless pair 
interaction. Quite arbitrarily we will say that the interaction is attractive if the 
Fourier transform w(q) is positive for all q; in the converse case it will be said 
repulsive. 

In a given GC configuration C = (N\ x\ . . . xn) the microscopic density of parti- 
cles reads 



N 

p(x\C) = Y / S {3 \x-x l 
and the GC partition function E [u] can thus be written as 
E\iA = Tr 



cxp ( -PVhs(C) + ~ (p\w\p) + (V\p) 



OO 1 

Tr[...] = £^7 M-..rfn... 



N=0 JV ' 



where i = x^ and di = d?Xi. For a given volume V, E [u] is a function of (3 and a 
convex functional of the local chemical potential v(x) which we have strength- 
ened by using a bracket. In equation (J2J), exp(— /3Vhs(C)) denotes the hard sphere 
contribution to the Boltzmann factor and V = v + vs where vs — —w(0)/2 is f3 times 
the self-energy of a particle. It follows from our assumptions on w(r) that v$ is a 
finite quantity which depends however on the regularization of the potential in the 
core. In the right hand side of equation (J2J we have introduced convenient brac-kets 
notations 

(V\p) = f dl V(l)p(l) (3) 
Jn 
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(p\w\p) = Jjld2p(l\C)w(l,2)p(2\C) ■ (4) 
We show now how to re-express H [v] as a functional integral. 

2.2. The Hubbard-Stratonovitch transform 

The Hubbard-Stratonovich transformation [ |H1 E] which was proposed in the 
50ies is based on simple properties of Gaussian integrals. In this framework the GC 
partition function S [u] can be formally re-expressed as a functional integral [ |3J HJ 
IT3*1 rjJl rT5J ri^]. Moreover H [i/], when re- written in this manner, can be interpreted 
as the GC partition function of a fluid of bare hard spheres in the presence of an 
external random Gaussian field tp with a covariance given by the pair potential [ 
r^rjTjriH]. It will be convenient to distinguish the case of attractive interactions 
(w(q) > 0) from that of repulsive ones (w(q) > 0). 

• i) Attractive interactions (w(q) > 0) 

We start from the well-known property of Gaussian integrals 

1_ \ JVip exp (-§ ((flw-^ip) + (p\(p) 



ex P (2 (P\ W \P) 



N w = JV(p exp (-± (Lp\w- l \(p)^ 
(aq>«pl<p))) w , (5) 



where ip is a real random field and the inverse w~ l must be ununderstood in 
the operator sense, i.e. J n d3 u>(l, 3)w~ 1 (3, 2) = 5(1,2). The functional inte- 
grals which enter the above equation can be given a precise meaning in the 
case where the domain Q is a cube of side L with periodic boundary condi- 
tions (PBC) which will be implicitly assumed henceforth. Additional technical 
details are to be found in the appendix. We now insert equation (J5J) in the 
definition (J2J of the GC partition function H [v] which yields 

~M = (Tr[exp(-(3V HS (C) + (V+<p\d)]) w , 

= (E m [u + p]) w . (6) 

S [u] , when re- written in this manner, can thus be re- interpreted as the GC 
partition function of a fluid of bare hard spheres in the presence of a random 
Gaussian field ip with a covariance given by the pair potential [ flT H IT7 | IT8] . 

ii) Repulsive interactions (w(q) < 0) 
This time we make use of 

/_1 \ /£V exp ((plw^lcp) + i(p\(p)) 

^ (exp(i(p|^)))^ , (7) 
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which yields 

SM = (S H s[I7 + ^]) ( _ w) • (8) 

Note that since (p is a real scalar field Shs [v + has to be evaluated for 
imaginary chemical potentials which may cause some troubles since log Shs H 
has singularities (i.e. branch cuts) in the complex plane [EH- In the repulsive 
case, the hard core part of the interaction is in fact not compulsory to ensure 
the existence of a thermodynamic limit [|2D1 an d the reference system can be 
chosen as the ideal gas [ HHJ |2I| • 

Some comments are in order. Firstly, equations (jUJ) and (jHJ) are easily generalized to 
the case of mixtures [^J|22| or molecular fluids. Secondly, when the pair interaction 
w is neither attractive nor repulsive, it is necessary to introduce two real scalar fields 
(f + and if some rigor is aimed at [ llOj . Alternatively, in a more sloppy way, eq. © 
can be considered to hold in any cases having in mind that <p will be a complex scalar 
field in the general case. We will thus write formally in all cases 

2 M = AC 1 / T><p exp {-H K [v, <P\) , (9) 
where the action (or effective Hamiltonian) of the KSSHE field theory reads as 

Ti-K [v, p] = ^ (^k _1 |y?) - ln5 HS P + (p] ■ (10) 

2.3. The collective variables representation 

We introduce now briefly the CV representation of H [u] and refer the reader to 
a vast literature for a more detailed presentation (see, e.g. [0 El 03 1221 121] )• The 
starting point is the formally trivial identity 

ex P Q (PIMP)) = J v p &f \p - P) ex P Q (pMp>) , (ii) 

where 5? [p] denotes the functional "delta" [|2]. Making use of its functional integral 
representation (see the appendix) 

6r[\] = Jt>u exp(z(o;|A)) , (12) 

one finds for the GC Boltzmann factor 

ex P Q (PHP)) = J VpVu exp (- (p\w\p) + i (u\ {p - p})^j . (13) 

Inserting equation (fT3|) in the definition (J2J) of the GC partition function 5 \y] one 
obtains 

E[u] = J Vp exp Q (pMp>) J\pA , (14) 

where the Jacobian 

J [p,u] = Vu exp (i (u\p)) S H s P - iu] (15) 
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allows for the passage from the microscopic variables %{ to the collective ones p. We 
note that J [p,v] does not depend on the pair interactions w(l, 2) but only on the 
GC partition function of the reference HS system S H s t* 7 ] which is supposed to be 
known. 

Equation (|14j) can also easily be recast under the form of a standard statistical 
field theory, i.e. as 

E[v] = J VpVuj exp (-Hcv [v, p, uj}) , (16) 

where the action of the CV field theory reads as 

Hcv[v,p,u] = -- {p\w\p) -i(u\p) -lnH HS [v-iu] . (17) 

We stress that uj and p are two real scalar fields and that equations (fTBj) and (JT7|) 
are valid for repulsive, attractive as well as arbitrary pair interactions. Moreover, 
with the clever normalization of Wegner for the functional measures there are 
no unspecified multiplicative constant involved in equation ([To]) (see the appendix 
for more details). 

The CV transformation is clearly more general than the KSSHE transformation 
since it can be used for a pair interaction w(l, 2) which does not possess an inverse 
and is easily generalized for n-body interactions (n > 2). The equivalence of the 
CV and KSSHE representations © and (fTB^l of H [u] is readily established in the 
repulsive case ( w(q) < 0) by making use of the properties of Gaussian integrals (cf. 
equation (J7|)). In the attractive or in the general case we cannot propose a convincing 
way (i.e. a non formal one) to establish this equivalence however. 

2.4. Statistical average 

In the sequel it will be important to distinguish carefully, besides the usual GC 
average < A(C) >gc of a dynamic variable A(C), between two types of statistical 
field averages. At first the KSSHE averages defined as 

(A M) K = 2 [u]- 1 JvpA [p] exp (-H K [u, <p]) , (18) 

where A [<p] is some functional of the KSSHE fields ip and, secondly the CV averages 
defined in a similar way as 

(A [p, uj]) cy = S [v]' 1 J VpVoj A [p, uj] exp {-H C v [v, P, w]) , (19) 

where A [p, uj] is an arbitrary functional of the two CV fields p an uj. 

3. Correlation functions 

To establish the link between the usual theory of liquids and the two statistical 
field theories which were introduced in sections (J2.2)) and (J2.3)) means to find the 
relations between density correlation functions in the one hand and field (either 
KSSHE or CV) correlation functions on the other hand. This is the purpose of the 
present section. 
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3.1. Density correlations 

The ordinary and truncated (or connected) density correlation functions of the 
fluid will be defined in this paper as [ [TTJ HZ] 

G (n) H(l,..,n) = (f[p(x t \C)\ 

\l=l / GC 

f 5 n S[v] 
E[u] 6u{l)...6u{n) ' 

G^ T M(l n) = ^ lQgSM (20) 

Our notation emphasizes the fact that the G^ n ' (connected and not connected) are 
functionals of the local chemical potential u(x) and ordinary functions of the space 
coordinates (1, . . . , n) = (x\, . . . , x n ). We know from the theory of liquids that [ 

EH dZ] 

GM' T M(1, • • • ,n) = GW[«/](1, . . . ,n) - E II G^^WKh, ■ ■ ■ ,hn) , (21) 

where the sum of products is carried out over all possible partitions of the set 
(l,...,n) into subsets of cardinal m < n. Of course p[^](a;) = G^ n=1 ^](a;) = 
G^ n=1 '' T [u](x) is the local density of the fluid. 

It follows from the definition ()20|) of the G ( - ri ^[z/](l, . . . , n) that they can be reex- 
pressed as KSSHE or CV statistical averages, i.e. 

G (n) M(l,..,«) = (GfeV+d(l,-..,n)) K , (22) 

G (fl) M(l.-,n) = (4 n sV-H(l,-,4v ^ 
Although enlightening these relations are not very useful except for the special case 
n — 1 which reads explicitly as 

p[u)(x) = (f>Hs[V+(p](x)) K , (24) 
p[v](x) = (pBs[?-iw](x))cv > ( 25 ) 

where Phs [£](£) is the local density of the hard sphere fluid at point x in the presence 

of the local chemical potential 

3.2. Field correlations 

The correlation functions of the KSSHE field ip and the CV fields p and uj will 
of course will be defined as 

G< n) H(l,..,n) = (flvixi)) , 

\l=l / K 

GM[u)(l,...,n) = (flp{xi)) , 

\l=l / cv 

G<WM(l,...,n) = (ft , (26) 

\i=i / cv 
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and their truncated counterparts will be defined in the same way as in equation ()21)1 . 
The relations between the field correlation functions ()26|) and the density correla- 
tion functions G^ n ^](l, . . . , n) are easily obtained by introducing ad hoc generating 
functionals and we quote here only the results, referring the reader to references [ 
ITUl 125] for a detailed dicussion. 

• i) relations between and G^ n > 
One has 

(<p(l)) K = w(l,l')p[u)(l'), 
Gf' T M(l,2) = w(l,2)+w(l,l')w(2,2')G^ T [»}(!', 2'), 
GW' T H(l,...,n) = w(l,l')...w(n,n')x 

x G (n)>T [v\ (1', . . . , n') for n > 3 , (27) 

where we have adopted Einstein's convention, i.e. space integration of variables 
labeled by the same dummy indice is meant. It will be quite enlighting to 
consider the GC correlations of the microscopic potential 

N 

V(x,C) = Yw{x,Xi) = / d 3 x p(x',C))w(x,x') . (28) 

We deduce readily from equations (f2Tjl that 

(^(1)) K = <^(D) GC , 
Gf'>](l,2) = (v(l)V(2))l c + w(l,2), 

GW' T [i/](l,...,n) = (V-(l)...y(n))^for n > 3 . (29) 

Therefore the KSSHE field (p(x) identifies "nearly" with the microscopic field 
V(x) at a neutral point. 

• ii) relations between G)?\ G^> and 
One finds easily that 

G^ T [u](l,...,n)=G^ T W](l,...,n), (30) 

which was expected and that 

G^ T [v] (1, . . . , n) = {-iY G^ T [v] (1, . . . , n) , (31) 

which was suspected. We stress that the above relations are valid for all n. In 
particular, for n — 1 note that (p(x)} CY = p(x) and (a>(x)) cv = i(y(x)^ Gc . 
The CV p and a> are thus connected to the density of particles and the density 
of energy respectively. We have therefore the correspondance ip <-> — iuo which 
will be a guideline in the sequel. 
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4. Mean Field theory 

Contrary to the situation which prevails in the theory of liquid, the term "mean- 
field" approximation is defined unambiguously in statistical field theory. It consists 
in neglecting field fluctuations and in approximating the grand potential (i.e. the 
logarithm of the partition function) by minus the action at its saddle point [Q] . In this 
section the KSSHE and CV theories will be considered at this level of approximation. 

4.1. KSSHE representation 

The mean-field (MF) or saddle point approximation of the KSSHE theory will 
therefore be defined by the set of relations [ E| : 

2 M f H = exp (-Hk [v, <Po]) , (32) 
where, for = (p , the action is stationary, i.e. 

S H K [u, if] 



5ip 



. (33) 



Replacing the KSSHE action by its expression (jlPji in eq. (J3~3*)) leads to an implicit 
equation for ip : 

<p Q {l) = w(l,l')p HS [V+<p ](l'), (34) 

which reduces to 

ip = w(0) pus P + <Po] (35) 

for a homogeneous system. It follows from the stationary condition ()33|) that the 
MF density is given by 

Pmf [v\ (1) = ^=MlM = PHS [v + (1) , (36) 

and that the MF grand potential reads 

lnS M F [v\ = lnH H s \?+ <^o] - ^ (pmf\w\pmf) ■ (37) 
Moreover, the MF Kohn-Scham free energy defined as the Legendre transform 

P-Amf [p] = sup { (p\ v) - In E MF [u] } (38) 

is found to be the convex envelop of 

/?Aif [p] = PA m [p] - - (p\w\p) + lJ Q dx w(0)p(x) . (39) 

It can be shown [ [TU] that ^4mf [p] constitutes a rigorous upper bound for the exact 
free energy A [p] if the interaction is attractive (w(q) > 0) and a lower bound in the 
converse case (w(q) < 0). 
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Finally, the pair correlation and vertex (or direct correlation) functions at the 
zero-loop order which are defined respectively as 

g&m(i,2) = f °:r , (40) 

o^(l) d)z/(2) 

are given by 

G®(1,2) = (1-«;^G^[I7+^]) _ \g^[I7 + ^ ](1,2) (42) 
c£(l,2) = -Gm 1 f (1 ) 2) = C hs (1,2)+ W (1,2). (43) 

We note that these relations coincide with the old random phase approximation 
(RPA) of the theory of liquids [EJ. Taking advantage of the arbitrariness of w(r) 
in the core on can chose a regularisation of the potential such that the pair corre- 
lation function g{r) = for r < a, then we recover the ORPA (optimized RPA) 
of the theory of liquids [ [TT] . To be complete let us point out that it follows from 
equation (|2*7jl that the free propagator of the KSSHE field theory is given by 

G^ 2) MF (1, 2) = (l-w* OH P + Vo] (1, 2)) ^ * w(l, 2) . (44) 



V, MFV 

4.2. CV representation 

An analysis similar to that of section (j4.1|) can be made in the CV representation. 
The MF level of the CV field theory will be defined by 

Smf [v] = exp {-Hey W, Po, wo]) , (45) 
where, for p = p and to = u the CV action is stationary, i.e. 



<$ ^cv [v, P, w] 



5p 



5 Hew [v, p, u] 



= . (46) 

(po,^o) 



Replacing the CV action by its expression (fTTj) in equation (|47)j) yields a set of two 
coupled implicit equations for p and u : 

= w(l,2)po(2)+iw (l) , 

= Po(1)-Phs[^-^ ](1) • (47) 

If we define </3 = — zu;o, then the two previous equations can be rewritten 

<p = w(l,2)p (2) , 

Po(l) = Phs^+<A>](1), (48) 

which shows that, as expected, </?o coincides with the saddle point of the KSSHE 
field theory (cf Sec 14. Moreover a direct calculation will show that 

In S M f [v] = -Hk [v, <fo\ = -H C v [v, p , u ] . (49) 

Therefore the local density, the grand potential, the free energy, the correlation and 
the direct correlation functions coincide at the MF level in the CV and KSSHE field 
theories. 
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5. Loop expansion 

The loop expansion is a usefull tool of statistical field theory [[I]. Away from the 
critical point it yields approximations for the free energy beyond the MF level. 

5.1. Loop expansion of the grand potential 

The loop expansion of the logarithm of the partition function of a scalar field 
theory can be found in any standard text book, see e.g. that of Zinn- Justin [[I]. In 
the case of the KSSHE field theory we have just to reproduce known results. One 
proceeds as follows. A small dimensionless parameter A is introduced and the loop 
expansion is obtained by the set of transformations 



In 2 \v 



(p = (p + A 1/2 x , 



AlnjAC 1 / V X exp 



A )) 1 

In E<°> [u] + A In S« [v] + A 2 In [ u ] + £>(A 3 ) , (50) 



where ipo is the saddle point of the KSSHE action. The expansion of In S [v] in powers 
of A is obtained by performing a cumulant expansion and by making a repeated use 
of Wick's theorem. At the end of the calculation one set A = 1. 

In order to obtain a loop expansion of InS in the CV representation a little 
homework is necessary. In this case we are led to introduce the following set of 
transformations 

p = po + X^Sp, (51) 
u = co + \ 1/2 5lu , (52) 

In ~ [„] = Am j/ V5p V5uj exp ^ev^El^ j J , (53) 

where (p ,u ) is the saddle point of the CV action. The form retained in equa- 
tions (|SH) and (JH2|) is imposed by the exact relation (u;} cv (1) = iw(l, 1 ) (p) cv (1 ) 
(cf. equations (j3T)j) and (}3*T]) for n — 1). Then one performs a cumulant expansion 
of (p)3*j) and use the appropriate version of Wick's theorem (i.e. for two coupled 
Gaussian fields p and u). The derivation is detailed in reference [|2E] and will not 
be reproduced here. One finds that the KSSHE and CV loop expansions of In S [u] 
coincide at each order of the expansion (i.e. the lnS^ [u\ are the same) which was 
expected. 

Of course S^ ) [u] = Smf [v\ ■ At the one- loop order one finds [ HUl P2H] 



s(1) [r] _ -A/"a, _ Jgg exp (-1(^1^ ^ (54) 
Nv> JVip exp (-|(y?H^ 



where we have adopted Zinn- Justin's notations for the free propagator [ri] 

A-Vl, 2) = w-\l, 2) - OH' T P + <Po] (1, 2) . (55) 
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Of course A ¥ , (l,2) coincides with G^ MF (1,2) as can be seen by comparing equa- 
tions (|44j) and (J55j) . If v is uniform then the system is homogeneous and A ¥ , (l,2) 
takes on a simple form in Fourier space, i.e. 



w{q) 

l-w(q)G^ T [V+ Vo ] (q) 



and, in this case, the Gaussian integrals in equation (p>4"j) can be performed explicitly 
(cf the appendix for more details). One finds that 

lnS«[i/] = -£ Jfo(l-w(q)G™' T [V + <p ](q)) . (57) 

The second-loop order contribution In S^ 2 ) to the grand potential has a complicated 
expression involving the sum of three diagrams sketched in figure. ^ 

lnS^M = D a + D b + D c . (58) 




D a D b D c 

Figure 1. Diagrams which contribute to lnH^ 2 ) [u], D a and D c are irreducible 
while D is reducible. 



More explicitly one has [[I] 

D a = I [dl...d4 A W) (1,2)A^ (3,4)GS T [V + <p ] (1,2,3,4) 

o JQ 

D b = I f dl...d3dl ...d3 A^ (l,2)A^ (l',2')A V0 (3,3') 

8 JQ 

x Gg' T [V + (1, 2, 3)Gg' T [77 + ^ ] (1', 2', 3') 
Dc = ^/ dl...d3dl'...d3' A w (l,l')A w (2.2')A Wl (3,3') 

x Gi' T [l7 + ^ ](l,2,3)Gi' T [l7 + ^ ](l',2',3'). (59) 

As they stand, the above relations are not particularly useful for practical applica- 
tions (even in the homogeneous case) since they involve the three and four body 
correlation functions of the reference HS fluid. We will introduce some reasonable 
approximation in section ()5.2)1 to tackle with this horrible expression. Quite remark- 
ably it has been shown recently that for a symmetric mixture of charged hard spheres 
lnH^ 2 ) has a much more simple expression which involves only the pair correlation 
functions of the HS fluid as a consequence of local charge neutrality [ HH1 1221 • 
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5.2. Loop expansion of the pressure and the free energy 

In this section we restrict ourselves to the homogeneous case, therefore In S [u] = 
VpP^u), where P denotes the pressure and PA [p] = V/3f(p) where / is the Hel- 
moltz free energy per unit volume. The two- loop expression that we derived for P 
in section (J5.1|) is too much complicated to be of any practical use since it involves 
the 3 and 4 body density correlation functions of the HS fluid which are unknown, 
whereas Gjjg is known approximatively, for instance in the Percus-Yevick (PY) 



approximation [ rn] . A simple but coherent approximations for G^g T and G 
will be proposed now. 

Recall first that it follows from their definitions [1213121] ( see e -S- equations (f2*U|) 
that the G^g' T [v] satisfy to the following relations 



(4),T 
r HS 



H« + i) Ghs 



(n),T 



V\ 



G 



(n+l), T 
HS 



n, n + 1) . 



(60) 



For a homogeneous system (in which case v is a constant) one infers from this 
equation that 



J dl...dn G^ s +1) ' T [u] (1, . . . , n, n + 1) = ^Phs W) 



(n) 

Phs 



~ (61) 

where pns {v) is the number density of hard spheres at the chemical potential v. In 
the rest of the section we will adopt the following approximation 



,71 



(n+l) 

Phs 



» S(n, 1) . . . 5(2, 1) for 72 > 3 



(62) 



Note that this hypothesis is coherent with the exact relations (JoTj|) and (J61|) . However, 



we keep the full G^g T [v\ (1, 2) = G^ s ' 1 [u] (xi 2 ) with the assumption that it is a 
known function of x±2 (in the PY approximation for instance). The free propagator 
A lfi0 (xi2) has thus also a well defined expression. It is not difficult to convince oneself 
that that the set of approximations that we have introcuced is reasonable as long as 
the range of the KSSHE field correlation functions is (much) larger than than the 
range of the HS density correlation functions. This will be true if w is a long range 
pair interaction (near the Kacs limit for instance). 

With the hypothesis it is easy to obtain the two-loop order approximation 
for the pressure. One finds 



t(2),T 



(3P {0 \u) 



pP {0 \u)+\pP 



(i)/ 



\ 2 f3P {2 \v) + C(A 3 ) 



PPhgpty) = Phs(v + <Po) - 



<Po 



2w(0) 



f3P 



\n(l~w(q)G^ T [V+ Vo ] (q) 



(3) 

^ o (0) + 



' (2) 

Po 



+ 



(2) 

Po 



2 



-^PoWKa (°) 



12 



d 3 x Aj ( 



x) 



(63) 
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where p^ ] = p^ (u + <p ) and J q = J d 3 q/(2n) 3 . 

The free energy (3f(p) is obtained as the Legendre transform of (3P(v). The 
calculation is lengthy and tricky and we must contend ourselves to give the final 
result (see reference [|2E| for the complete proof). Setting A = 1 one finds 



Pf(p) 



PfBs(j>) T^P 



+ - 



In (l-w (q) G 



1 

P-A (0) + -A (0) 



T 

HS, p 
7 (2)' 



(?) 

2 



pw (?)} 

n(2) n2 



)(!) 



12 



d 3 x Aj(x) 



(64) 



where the HS pair correlation function G^ s must be evaluated at the density p 



and, similarly, 

Finally, in equation (fBIjl 

n( 2 ) - 



Ap(<?) 



g hs, P (?) ™ (g) 



(65) 



(i) 

HS 



^ _ -»miP) 



(2), 



9^ 



fflw 



p 



(3) 





3 


^I(P)" 


2 ~ I/ ffi(p) I/ Hs(p) 


dv 






"41 (P)" 


5 



(66) 



where ^hs(p) denotes the nth derivative of the HS chemical potential with respect 
to the density (it can be computed in the framework of the PY or Carnahan-Starling 
approximations for instance [E])- Some remarks are in order. 

i) It must be pointed out that, quite unexpectedly, the reducible diagram Df, 
has not be cancelled by the Legendre transform. Usually, in statistical field theory 
it is the case (cf. [HJ). The reason is that the chemical potential v is not the field 
conjugate to the order parameter m = (<p) K of the KSSHE field theory. However 
one of us have shown elsewhere [ 122 that for the symmetric mixtures of charged 
hard spheres only irreducible diagrams contribute to (3f^\ 

ii) The two first terms of equation (|6"4"j) (i.e. the one-loop order result) supple- 
mented by equation (|44j) for the pair correlation function are exactly constitutive of 
the RPA theory of liquid. The full result including the two-loop contribution is a 
new approximation as far as the authors know. 

iii) All the quantities which enter equation ()64j) can be computed numerically 
(for instance in the PY approximation). However we stress once again that the 
validity of equation is restricted to long range pair potentials (for instance 
w(r) = 7 3 $(r/7) with 7 — > 0) and we are not aware of exact results (i.e. obtained by 
Monte Carlo simulations) which would allow us to test the validity of our expression 
for the free energy for such kind of fluids. 
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6. Conclusion 

Using the CV method we have reconsidered the basic relations of the statistical 
field theory of simple fluids that follow from this approach. In contrary to the KSSHE 
theory [ Elj the corresponding CV action depends on two scalar fields, the field 
p connected to the number density of particles and the field u conjugate to p, 
connected to the density of energy. Explicit relations between field and density 
correlation functions were obtained. 

For a one-component continuous model of fluid, consisting of hard spheres inter- 
acting through an attractive pair potential, we have calculated the grand partition 
function for both versions of statistical field theory using the loop expansion tech- 
nique. As it was expected, at all the orders of loop expansion considered, both 
versions of the theory produced indeed the same analytical results. The expressions 
for the pressure as well as for the free energy were derived at the two-loop level. It 
yields a new type of approximation and Monte Carlo simulations are wanted to test 
its validity. 

From our analysis of the CV and KSSHE transformations we can also conclude 
that the former has some important advantages which could be very useful for more 
complicated models of fluids. In particular, it is valid for an arbitrary pair potential 
(including a pair interaction w(l, 2) which does not possess an inverse) and is easily 
generalized for the case of n-body interparticle interactions with n > 2. 

The statistical field theory of liquids is still in its enfancy and many progresses 
including the implementation of Wilson's ideas on the renormalization group for 
studying gas-liquid critical points are likely to be achieved in the future, this is at 
least the hope of the authors. 
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Functional measures and integrals 

We give in this appendix some details concerning functional measures and in- 
tegrals. Let us consider a real scalar field <f(x) defined in a cube C 3 of side L and 
volume V = L 3 . We assume periodic boundary conditions, i.e. we restrict ourselves 
to fields which can be written Fourier series, 

where A = (2tv/L) 1? is the reciprocal cubic lattice. The reality of ip implies that, 
for q ^ <p q = (p*L q , where the star means complex conjugation. Following Wegner [ 
12*5] we define the normalized functional measure T>(p as 

Vcp = TT Jh= (68) 
d(p q d(f>-q = 2 dRefiq dlmipq for q ^ . (69) 
Eq. ([""{ can thus be conveniently rewritten as 

= TT dRe ^ dlm ^ , (70) 

where the sum in the right hand side runs over only the half A* of all the vectors of 
the reciprocal lattice A (for instance those with q x > 0). The normalization constant 
of equation (JHJ) can thus be explicitely evaluated as 



Af w = J Vip exp (-- (ip\w V) 

= ex P [\ H lnt %)J ~ ex P / ln *%)) » ( 71 ) 



where 2) is positive and satisfies w(l,2) = w(2, 1) = w(xi 2 ) and (27r) 3 / (? = 
/ (i 3 g. It is worth noting that in this case we have trivially Af w -i = 1/Af w . More 
generally and with the same hypothesis we have the two useful identities 

(exp ((z) (<p\u))) w = AC 1 / V<p exp (-± (p|«rV) + («) (</^>) 

= exp(V(-)^(u;|H^)) , (72) 

where a; is a real scalar field. 

We define now the "functional delta" distribution 5? [A] as 

6jr[X] = jvoo exp(i(u;|A}) , (73) 
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where both to and A are real scalar fields defined on C. Since 

HA) = f d 3 x cu(x)X(x) = ^J2 C -A > ( 74 ) 
^ V q eA 

it follows from equation that we have more explicitly 

5jr [A] = V2ttV 5 (a ) J] [ttV" <J (ReA,) 8 (lm\ q )] . (75) 

<?eA* 

Therefore 

Jv\5 F [\] = l, (76) 

and, more generally 

|l?A^[A]^[A-Ao]=^[Ao] , (77) 
where T [A] is some arbitrary functional of the field A (a;). 
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